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Abstract 

We propose an interpolation method that is invariant under 
Mobius transformations; that is, interpolation followed by 
transformation gives the same result as transformation fol- 
lowed by interpolation. The method uses natural (Delau- 
nay) neighbors, but weights neighbors according to angles 
formed by Delaunay circles. 

1 Introduction 

An inversion of R 2 U {00} is a bijection generated by a 
circle, that maps circles to circles and preserves angles 
between pairs of curves. The circle with center c and ra- 
dius r gives an inversion as follows: each point p along a 
ray with origin c is mapped to another point p' along the 
same ray, so that the product of the distances \cp\ and \cp'\ 
is r 2 . Thus points on the circle remain fixed and the cen- 
ter c maps to 00. The set of products of inversions forms a 
group, the group of Mobius transformations on Euclidean 
space E 2 . 

Mobius transformations arise in data visualization. For 
example, hyperbolic browsing [^j, || uses the Poincare 
model of the hyperbolic plane to view large graph struc- 
tures. A change of focus in a hyperbolic browser is 
given by an isometry of hyperbolic space correspond- 
ing to a Mobius transformation of E 2 that maps the unit 
disk to itself Brain mapping [^j flattens a convo- 
luted surface by mapping it quasiconformally to a disk; 
subsequent navigation in the brain map can be accom- 
plished by hyperbolic browsing. For these applications, 
it is reasonable to ask for computational primitives spe- 
cially designed for compatibility with Mobius transforma- 
tions. We previously considered the problem of finding 
an optimal Mobius transformation for viewing hyperbolic 
or spherical data [[j]]. In this paper we give a Mobius- 
invariant method of extending a function defined only at 
discrete sample points to a continuous function defined 
everywhere on their convex hull. 
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2 Natural Neighbors 

Let S = {si,S2, ■ ■ ■ , s n } be a set of points in E 2 , and let s 
be another point in the convex hull of S. Assume we have 
real- or complex-valued "elevation" zi for each s,, and we 
would like to define an elevation z =f(s) at s, such that/ 
is a continuous function of s and/(s,) = z,-. 

Natural neighbor interpolation [Q| sets/(s) to be a con- 
vex combination of the Zi values at the natural neighbors 
(also known as Delaunay or Voronol neighbors) of s. Let 
V denote the Voronoi cell of s in the Voronoi diagram of 
S U {*}) (that is, the set of points in E 2 closer to s than 
to any si), and let V, denote the Voronoi cell of St in the 
Voronoi diagram of S. Then/(s) = J2i w iZi where w; is 
the fraction of V covered by V,-. See Figure |l](a). Sib- 
son [Q] showed that the weights w, act as a sort of local 
coordinate system for s in the sense that if s = (x,y) and 
■s; = (xi,yi), then Xi = £\ w,*,- andy,- = J^w,-?,-. Sibson's 
theorem also implies that natural neighbor interpolation 
reconstructs linear functions, that is, if z; = axi + by, + c 
for each i and some a, b, and c, then z = ax + by + c. 

Natural neighbor interpolation is a logical starting point 
for Mobius-invariant interpolation, because the notion of 
Voronoi neighbor can be made invariant under Mobius 
transformations. Two sites Sj and sj are extended Voronoi' 
neighbors if there is a circle C through s, and sj bounding 
either an empty disk or empty disk complement. After an 
inversion, C or its complement remains empty and thus s, 
and sj are still neighbors. 

3 What Weights? 

What do we do with the weights wf! Area changes un- 
der Mobius transformations, so we must find an invariant 
alternative: angle. The angle most obviously associated 
with natural neighbor s, is the angle at the lune formed 
by maximal empty circles, as shown in Figure [jjb). We 
use the exterior lune angles 6*,-, because the sum of these 
angles is fixed: Yli ®i = 27r. (One way to see this fact 
is to perform an inversion that takes s to 00; the exterior 
lune angles become the exterior angles around the convex 
hull of the transformed point set.) For invariance, we can 
simply set w, to be proportional to any function of angle, 
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Figure 1: (a) Classical natural neighbor interpolation weights neighbors according to fractions of the area of the 
Voronoi cell of s. (b) Mobius-invariant natural neighbor interpolation weights neighbors according to lune angles. 



w(9j). But is there some especially natural choice for w()? 

First, we would like the interpolation to be continuous 
as s — ► Si, so we require w(8j) — > oo as 6*,- — * ir. Sec- 
ond, we would like the interpolation to reconstruct some 
class of functions, analogous to linear functions in Sib- 
son's method. Harmonic functions (solutions to Laplace's 
equation V 2 f = 0) would be high on anyone's wish 
list, because of numerous applications (heat, electrostat- 
ics) and because the family of harmonic functions is itself 
invariant under Mobius transformations. The values of a 
harmonic function on a discrete point set, however, do not 
determine the function, so we cannot hope for a property 
as strong as Sibson's theorem. We instead ask for a prop- 
erty that holds only in the limit. 

Lemma 1. Let s\, S2, ■ ■ ■ , s n be a set of points on a circle 
C, and let e be the maximum angle of arc between a pair 
of successive Si's. Let s be a point interior to C, and let 
9i be the exterior lune angle between s and Si as above. 
Now let f be a harmonic function defined on the closed 
disk bounded by C. Then as e — * 0, ^ ~~ > f( s )- 

Proof sketch: There is a Mobius transformation that puts 
the image of s at the center of the image of circle C. Now 
the harmonic measure of an arc of a disk with respect to 
the disk center is proportional to its arc angle (see for ex- 
ample [jfij]). It is not hard to confirm that in the limit of 
small arcs, the arc angles are the same as the exterior lune 
angles. Thus J2i ®if( s i) i s the sum over small arcs of har- 
monic measure times the value of/() at a point in the arc; 
this sum converges to the limit /(*). ■ 

Thus in summary we would like the weighting func- 
tion w(0j) to go to infinity as 0,- — > ir and w(9i)/9i to 
go to a constant as 9i — > 0. The most obvious choice 
is to make w(9i) proportional to tan(0,-/2). We normal- 
ize the weights so that they sum to one, and thus the 
actual weights are wt = (1/W) tan(0,-/2), where W — 
Eitan(0 ( -/2). 



Theorem 1. Natural neighbor interpolation using 
weights proportional to tan(0,-/2) gives a continuous 
function f(s) that interpolates the elevations at sample 
points and, in the limit of dense samples on a circle, 
reconstructs harmonic functions. 

Open questions include: Is there a nice C 1 continuous 
version? (One idea is to weight by tan 2 (#,/2), but this 
loses the harmonic function property.) What about higher- 
dimensional Mobius-invariant interpolation? 
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